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Introduction: Size distributions of nebular solids 
in chondrites suggest an efficient sorting of these early 
forming objects within the protoplanetary disk. The 
effect of this sorting has been documented by investi- 
gations of modal abundances of CAIs (e.g., [1-4]) and 
chondrules (e.g., [5-8]). Evidence for aerodynamic 
sorting in the disk is largely qualitative, and needs to 
be carefully assessed. It may be a way of concentrating 
these materials into planetesimal-mass clumps, perhaps 
100’s of ka after they formed. A key parameter is 
size/density distributions of particles (i.e., chondrules, 
CAIs, and metal grains), and in particular, whether the 
radius-density product (rxp) is a better metric for de- 
fining the distribution than r alone [9]. There is no con- 
sensus between r versus rxp based models. Here we 
report our initial tests and preliminary results, which 
when expanded will be used to test the accuracy of 
current dynamical disk models. 

Samples: Allende, a CV3 oxidized carbonaceous 
chondrite contains nebular components, including a 
diversity of chondrules and refractory inclusions (Fig. 
1A,B). NIST 1019b is a well-characterized particle 
size standard [10] that contains glass beads with di- 
ameters ranging from -750 to 2450 pm (Fig. 1C,D). 
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Figure 1. A) X-ray chemical map of Allende, with Mg, Ca, 
and A1 as red, green, and blue. B) Best fit ellipses defined by 
ImageJ for just porphyritic chondrules from A. C) An exam- 
ple radial section of the NIST 1019b particle size standard 
imaged by x-ray tomography [10]. D) Best fit ellipses de- 
fined by ImageJ of glass beads from C. 

Methods: The Allende meteorite was cut into -cm 
thick slabs and broken into fragments that could be 
made into one inch epoxy rounds. One fragment was 


sawed into four serial sections spaced -1.2 mm apart. 
The sections were polished using a series of grit sizes 
down to 1 pm for SEM analysis. 

We took x-ray maps of the Allende sample using a 
JEOL 7600 Field Emission SEM with a silicon drift 
detector. Images were taken at 15 kV, 30 nA, with a 90 
pm aperture. The maps were taken at ~150x magnifi- 
cation in which each pixel represents -3 pm These 
maps were stitched together to form a large montage of 
Allende to quantify the abundance and composition of 
the different nebular solids. Mg, Ca, and A1 x-ray maps 
were combined, and the particle distribution analyzed, 
using the image mixing and particle counting functions 
of ImageJ software. Initially, objects measured within 
the x-ray maps have been grouped textrually (i.e., po- 
phyritic, barred, radiating, etc.), but will ultimately be 
grouped by composition. Objects were measured for 
area using best fit ellipses. For simplicity we consider 
that the particles can be approximated by spheres and 
calculate an average sphere diameter from the major 
and minor ellipses axes measured for each. Diameters 
calculated in this manner match those determined from 
the circular area defined by the total number of pixels. 

We counted particles (i.e., nebular components) in 
Allende by tracing each using Adobe Illustrator. Then 
by taking the population of outlines and running it 
through ImageJ’ s automatic particle counting. We ob- 
tained the area, and major and minor axes of each par- 
ticle greater than -20 pm in diameter (seen in Fig. 2). 
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Figure 2. 2D size distribution of chondrules from -15% of 
the x-ray maps collected for the Allende CV3 chondrite. 

There are well-known sampling effects which 
cause the observed property N A (d), the number of par- 
ticle sections or profiles, of diameter d per unit area, to 
differ from the more fundamental, desired quantity 
N V (D), the number of spheres of diameter D per unit 



volume (e.g., [11]). Specifically, sections tend to cut 
particles non-diametrically, diminishing the fraction at 
true diameter D and artificially increasing the fraction 
at all smaller diameters (Fig. 3). We have written an 
algorithm that works as a simple matrix inversion [12] 
designed to unfold the particle sections obtained from 
the 2D Allende data set, however to calibrate the algo- 
rithm it was first tested on the particle size standard 
NIST 1019b. This was done using x-ray microtomo- 
graphy data [11] from 2D tomographic slices through 
NIST 1019b, a cylinder of glass beads. 

We counted the glass bead sizes in NIST 1019b in 
a manner similar to that for Allende. In our trials, glass 
beads were both hand counted, by drawing an ellipse 
around the bead and measuring that ellipse’s area, ma- 
jor, and minor axes, and auto counted by Image J, re- 
spectively (Fig. 1C,D). The error of measurements 
taken by hand were determined empirically to be 
-10%, which is similar to the uncertainty found by the 
several ways we determined diameter with the auto- 
mated Image J data collection routine. 



Figure 3. Graphical representation of the unfolding problem. 
The size distribution N A (d) and N V (D) are binned the same 
way, into a set of diameter bins with lower and upper 
boundaries at (di_ 1? di), with i = 1 representing the top end of 
the smallest particle size bin. For simplicity, we assume the 
spacing is regular, so di = di_i +A, with A = D max /N. 

Results and Discussion: The NIST 1019b standard 
was sampled near both the “top” and the “bottom” 
(sampling was made at intervals of 19.5 and 195 pm). 
The experiment was designed to capture the particle 
distribution of the overall population and differential 
settling within the sample container. The actual distri- 
bution and sampled distributions are shown in Fig. 4 A. 
Fig. 4B shows the data after it was processed through 
our unfolding algorithm. Our sampling does not re- 
solve the narrow peak at -800 pm. The peak is about 
50 pm wide, less than the absolute measurement error 
in the diameter estimate at that size. The success of our 
unfolding algorthm can be seen by comparison of Fig. 
4A,B from diameters of 0 to <800 mm where, after 
unfolding, all of the small diameters that reflect non- 
diametric slices through a large sphere have been cor- 
rected for (i.e., “Forbidden” sizes). The unfolded data 


(20-bin histogram) also resolves differential settling of 
smaller and larger size components. The first is shown 
clearly in the “top” subset of data at -1000 pm and the 
second in the “bottom” subset of data at -1200 pm. 

Allende x-ray maps revealed distinct chondrules, 
CAIs, and matrix in great detail, but future image 
processing will be needed to quantatively resolve more 
subtle intensity (i.e., compositional) differences. Nev- 
ertheless, preliminary analysis of -15% of the Allende 
data, shown in Fig. 2, demonstrate the general expo- 
nential decline of larger sizes (e.g., [13]). More data 
will be processed before meaningful density differ- 
ences among nebular components can be made. 
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Figure 4. Top panel shows the line histogram data generated 
from measurements of NIST 1019b. Bottom slice data is in 
blue, top slice data is in red, total slice data is in green, and 
the accepted values are in black. Bottom panel shows data 
after unfolding. Note that the settling shows up clearly as 
does the modes of the larger particles. 
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